Substantial reduction of Stone- Wales activation barrier in fullerene 
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Stone- Wales transformation is a key mechanism responsible for growth, transformation, and fusion 
in fullerene, carbon nanotube and other carbon nanostructures. These topological defects also 
substantially alter the physical and chemical properties of the carbon nanostructures. However, this 
transformation is thermodynamically limited by very high activation energy (~ 7 eV in fullerene), 
which can get reduced due to the presence of hydrogen, extra carbon atom, or due to endohedral 
metal doping. Using first-principles density functional calculations, we show that the substitutional 
boron doping substantially reduces the Stone- Wales activation barrier (from ~ 7 eV to 2.54 eV). This 
reduction is the largest in magnitude among all the mechanisms of barrier reduction reported till 
date. Analysis of bonding charge density and phonon frequencies suggests that the bond weakening 
at and around the active Stone- Wales site in B-heterofullerene is responsible for such reduction. The 
formation of Stone- Wales defect is, therefore, promoted in such heterofullerenes, and is expected 
to affect their proposed H2 storage properties. Such substitutional doping can also modify the 
Stone- Wales activation barrier in carbon nanotube and graphene naonostructures, and catalyze 
isomerization, fusion, and nanowelding processes. 

PACS numbers: 61.48.-c, 81.05.ub, 34.10,+x 



I. INTRODUCTION 

Stone- Wales (SW) defects^ are important topological 
defects in sp 2 -bonded carbon materials, and play cru- 
cial role in growth, isomerization, and nanoscale plastic- 
ity of carbon nanotube, fullerene, and graphitic nanos- 
tmcturesP^SW transformation involves an in-plane 90° 
bond rotation with respect to the bond center. This leads 
to pentagon-heptagon defects (5-7-7-5 dislocation dipole) 
in carbon nanotube and graphenep"^ and interchange of 
a pair of pentagonal and hexagonal rings in the fullerene. 1 
Such structural transformations are believed to be the 
fundamental processes for the coalescence of fullerene 6 
and nanotubePthe formation of molecular junctions for 
nanoelectronic devices,^ and plastic deformation. 3 Both 
chemically modified and unmodified SW defects induce 
local curvature to otherwise planer graphitic materi- 
als pKU w hich may enhance the formation of nanotube 
and fullerene from planer carbon nanostructures Co- 
alescence of fullerene and nanotube has been also pro- 
posed to occur through a sequence of such structural 
transformation.^ SW defects alter the electronic prop- 
erties of carbon nanostructures, and thus substantially 
modify chemical reactivity toward adsorbates (reactiv- 
ity increases compar ed to pristine counterpart),^"^ and 
transport properties ! 18 * 191 Similar SW defects are also 
observed in boron nitride nanotube and nanosheetj^ni 
and found to have important implications in determin- 
ing physical and chemical properties.^ 

It is known that 94% of the fullerene Ceo isomers can 
be derived by a sequence of SW transitions! 2 ^ Thus, 
SW transformation is believed to be the mechanism for 



fullerene isomerization or ih-Ceo formation, which rep- 
resents the global minimum on the pot ential energy sur- 
face, during high temperature annealing J 2 ^ 2 -^ Interchange 
of a pair of pentagonal and hexagonal rings in 2h-C6o 
through single SW transition leads to C2 V symmetry, 
which does not obey the isolated-pentagon rule. 25 This is 
the first isomer and separated by ~ 1.6 eV energy from 
the global 1^- minimum pSEH Therefore, C2V-C60 repre- 
sents the first step toward fullerene isomerization, or the 
last step of annealing process before reaching the icosa- 
hedral global minimum.^ Although the pristine C2V iso- 
mer is not experimentally observed earlier, the chlori- 
nated isomer has recently been stabilized experimen- 
tally via Kratschmer-HufTman synthesis} 2 ^ and in the gas 
phase by subsequent dechlorination.^ The SW transition 
in fullerene is, however, thermodynamically limite d due 
to very large activation barrier (~ 7 eV, Fig. [IJ)P 7 I 31 I 32 I 
The barrier height is found to get reduc ed in the pres- 
ence of an extra carbon or hydrogen atom! 2 -^ The extra 
carbon atom acts as autocatalyst and reduces activation 
energy by ~ 2 eV, when placed preferentially in the re- 
gions of paired pentagons. In contrast, the endohedral 
metal doping (K, Ca, and La) has found to be relatively 
less effective in barrier height reduction.^ 

Here we report the effect of substitutional doping on 
SW transition in Ceo through first-principles density 
functional calculations. We find that the activation bar- 
rier is reduced substantially by substitutional boron dop- 
ing at the active SW sites. The magnitude of reduction 
(2.12 eV and 4.34 eV for single and double B-doping, re- 
spectively) is much larger compared to the cases w ith ex- 
tra carbon or hydrogen or endohedral metal dopingP 27 * 33 * 
Boron doped heterofullerenes (Cqo- x B x ^ x ^ 6) have 
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FIG. 1. (color online) Activation energy barrier E a for /h-Ceo 
—t C2V-C60 transition, which represents single SW transition, 
is very high. Activation barrier for the reverse transition (CW- 
C60 — > i"h-C6o) is smaller by ~ 1.6 eV, which is the energy 
difference between Ceo in Ih- and CW- symmetry. The pyra- 
cyclene region is highlighted (red). 



been successfully synthesized} 3 -^ and are be lieved to be 
a promising candidate for hydrogen storage P^ 3 -^ It has 
been shown that H2 adsorption energy increases sub- 
stantially for B-heterofellerenes, where the B-centers act 
as adsorption centers. Our study shows that formation 
of SW defects are promoted in such B-doped hetero- 
fullerenes, which may alter the reversible storage proper- 
ties as SW defects are known to be more reactive towards 
adsorbates. Presence of such substitutional dopants is 
also expected to reduce the activation barrier and pro- 
mote SW transition in other carbon nanostructures, and 
thus catalyze isomerization, fusion, and nanowelding pro- 
cesses. 



II. METHODOLOGY 

Calculations were carried out using density functional 
theory implemented in the Vienna Ab-initio Simula- 
tion Package (VASP) code 37 with the Perdew-Burke- 
Ernzerhof exchange- correlation functional and projec- 
tor augmented wave pseudopotentiaP^. The kinetic en- 
ergy cutoff was chosen to be 800 eV. Symmetry unre- 
stricted geometry optimizations were terminated when 
the force on each atom was less than 0.005 eV/A, and 
the reciprocal space integrations were carried out at the T 
point. We determined the minimum-energy path for SW 
transition and the corresponding migration energy bar- 
rier E a using climbing-image nudged elastic band (NEB) 
method!^ In NEB, a set of intermediate states (images) 
are distributed along the reaction path connecting opti- 
mized initial (ih-C6o) an d final (C2V-C60) states. To en- 
sure the continuity of the reaction path, the images are 
coupled with elastic forces and each intermediate state 
is fully relaxed in the hyperspace perpendicular to the 
reaction coordinate. 



III. RESULTS AND DISCUSSIONS 

Determination of the activation energy necessary for 
SW transition (2h-C6o ^ CW-Ceo) via climbing-image 
NEB method requires structural information of the initial 
and final structures, within the concerned level of theory. 
Therefore, before we discuss the effect of B-doping on the 
activation energy, we begin our discussion with the struc- 
ture and electronic properties of Ceo, C59B, and C58B2 
cages. 



A. Ceo? C59B, and C58B2: Structure and Electronic 
properties 

There are two types of bonds in I^-Cqq: [6,6]-bonds 
(1.399 A) at the junctions of two hexagons are smaller 
than [5,6]-bonds (1.453 A) at the junctions of a pentagon 
and a hexagon. The optimized geometries in both 1^ and 
C2v symmetries are shown in Fig. [I] and 14-atom pyracy- 
clene region containing the rotating C2 dimer is shown in 
Fig. [2] Calculated [5,6] and [6,6] bond lengths (Fig.|2| are 
in e xcellen t agreement with previous theoretical calcula- 
tion d 32 * 4 ^ and experimental measurements^^ f or the ih 
structure. The rotating bond in C2V-C60 is much smaller, 
1.359 A, which is also in a greement with previous theo- 
retical calculations! 28 * 32 ^ However, no experimental in- 
formation is available on pristine C2V-C60, chlorinated 
C2V structure (CeoCls) shows a similar bond shortening 
of the rotating dimer (1.37 A)P3 As shown in Fig. [2] such 
bond shortening introduces local strain, and thus a dis- 
tribution of bond lengths is observed in the pyracyclene 
region. 

When one of the C atoms in the rotating dimer (see 
supplementary documentP^) is replaced by B, three bond 
lengths involving C and B are increased (two [5,6]- 
C-B=1.549 A, and [6,6]-C-B=1.524 A) for I h -C 59 B. 
This is responsible for the bond deformation in pyracyc- 
lene region (the [5,6]-C— C bond distribution ranges from 




FIG. 2. (color online) Local pyracyclene region for the op- 
timized Ceo in Ih and symmetries. The rotating dimer 
corresponding to the SW transition is highlighted, and the 
bond lengths are shown in A. The rotating bond length de- 
creases substantially in the Civ structure, and induces a sig- 
nificant local distortion. Similar trend is also observed in 
B-het erof ullerenes 
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FIG. 3. (color online) Orbital degeneracy for pristine and substitutional B-doped fullerenes (C59B and C58B2) for both ih and 
CW symmetries. The HOMO (LUMO) of the icosahedral Ceo, H u (Tin) is five (three) fold degenerate. Bottom (top) of the 
vertical arrow represents HOMO (LUMO) for the corresponding structure. Introduction of defects, both in terms of topological 
defect or substitutional B, breaks the orbital degeneracy, and thus reduces the energy gap E g (see Table III. 



1.442 A to 1.482 A). The trend of the structural changes 
in moving from ih-CsgB to C2V-C59B, remains same as 
is observed for pristine Ceo- Although, the structural 
properties of C58B2 show similar qualitative behavior, 
the pyracyclene region is comparatively more deformed 
due to the presence of long B— B (1.663 A) bond, and four 
longer C— B (1.551 A) bonds. Though, strictly speak- 
ing, both the icosahedral and symmetries are broken 
due to B-doping, for simplicity we shall retain the same 
nomenclature in the following. 

The structure is the first iso mer, a nd in agreement 
with previous theoretical reports) ^ 28 * 32 * we find that this 
structure is 1.57 eV higher in energy compared to the 
global ih-minimum. Calculated energy gap between the 
highest occupied molecular level (HOMO) and the low- 
est unoccupied molecular level (LUMO) is 1.63 eV for 
/h-symmetry, which agrees well with the spectroscopic 
measurements.^^ Due to the icosahedral symmetry, the 
HOMO (H u ) and LUMO (T lu ) are five and three fold de- 
generate, respectively, in Ih-C^o. These degeneracies are 
lifted (Fig. 3]) due to the presence of topological defect 
introduced by SW transition or due to the substitutional 
B-doping. Substitutional B introduces a local bond de- 
formation in the structure, and also introduces holes into 
the system. Both causes reduction of the HOMO-LUMO 



TABLE I. Energy difference between the ih and CW symme- 
tries, AE = E(Ih) — E(C 2w)\ and the energy gap E g between 
the highest occupied molecular orbital (HOMO) and lowest 
unoccupied molecular orbital (LUMO) for both ih- and CW- 
symmetry. All the energies are in eV. 





AE 


E 9 (h) 


Eg(C 2 v) 


Pristine-Ceo 


1.57 


1.63 


0.91 


C59B 


1.67 


1.18 


1.04 


C58B2 


1.75 


0.98 


0.70 



gap, which are shown in Table [T[ For pristine C2 V -^60^ 
the degeneracy of both HOMO and LUMO is completely 
broken (Fig. |3j), which in turn reduces E g by ~ 0.7 eV. 

Single substitutional B introduces a significant bond 
deformation,^ and thus the inherent symmetry of C59B 
deviates from perfect ih-symmetry, which removes the 
five-fold degeneracy of the HOMO level (Fig. [3). Sub- 
stitutional B introduces hole (one/B-atom) into the sys- 
tem, and thus the HOMO in I^-C^^B is singly occupied, 
which is pushed 0.47 eV above compared to the pristine 
Iyi-Cqo, which reduces the gap. In addition, when a SW 
topological defect is introduced the three-fold degener- 
acy of LUMO is completely lifted (Fig. |3|, and pulled 
down by ~ 0.2 eV. This reduces the gap even more for 
the C2V-C59B structure. Similarly, when two C atoms 
are replaced (the system has two holes) the degeneracy 
in HOMO is broken. One of these levels is completely un- 
occupied and serves as LUMO in ih-C58B 2 . This level lies 
0.98 eV above the occupied HOMO. The orbital degen- 
eracy of the HOMO is further broken and pushed above 
to further reduce the gap in C2V-C58B2 structure. 



B. Activation barrier: Effect of substitutional 
B-doping 

Activation barriers for the SW transition in carbon 
nanotube, graphitic nanostructure, and fullerene are gen- 
erally very high! 28 " 9 ! 50 ! Our calculated energy barrier for 
the pristine Ceo is found to be 6.88 eV (Fig. |4|, which 
is in agreement with previo us reports obtained within 
a range of theories! 2 ! 27 ! 28 ^ 32 ! Compared to carbon nan- 
otube or graphitic nanostructures, the activation barrier 
in fullerene is considerably less due to the strain present 
in fullerene structure. 28 We find that the rotating C2 
unit in the transition state (TS) pops-out (by ~ 0.21 
A) from the C60 surface. Similar TS was predicted pre- 
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FIG. 4. (color online) The activation energy barrier to the 
SW transition (lh CW) reduces substantially for the B- 
heterofullerenes, when B is doped in the active SW site (ro- 
tating dimer). The reverse barrier is smaller by an amount 
AE, reported in Table [I] which are similar in magnitude for 
all compositions. 



viouslyi ^ 1 32 1 51 1 The bond length of the rotating C2 unit 
shrinks to 1.248 A in the TS from the corresponding 
bond length in J h (1.399 A) and C 2v (1.359 A) symme- 
tries. Therefore, this C2 bond can be characterized as a 
stretched sp 3 -bondP^ 

We find that the large SW activation barrier is reduced 
substantially by B-doping (Fig. [4| at the active SW site. 
The barrier is reduced by 2.12 eV, when a single C-atom 
is replaced in the rotating dimer by B for the C59B het- 
erofullerene. The reduction is comparable in magnitude 
with the cases where an extra carbon was strategically 
placed in the regions of paired pentagons^ or catalyzed 
by hydrogen PS In contrast, the endohedral metal doping 
has relatively smaller effect. The barrier was found to re- 
duce by only 0.80 eV for La-doping In the present case, 
the barrier is further reduced to 2.54 eV, when both the 
C-atoms in the rotating dimer are replaced with B-atoms 
for C58B2 heterofullerene. Compared to pristine Ceo, a 
total reduction of 4.34 eV is obtained which is the maxi- 
mum reported till date. We find that for pristine Ceo as 
well as B doped Ceo, the structure corresponding to the 
saddle point has only one vibrational mode with imagi- 
nary frequency. This confirms that the obtained saddle 
points are indeed true first-order transition states. We 
find that the frequency of this imaginary mode decreases 
with the introduction of single B, which further decreases 
with the introduction of second B (962 cm -1 , 377 cm -1 , 
and 167 cm -1 for Ceo, C59B, and C58B2, respectively). 
Thus the variation of energy along the reaction path be- 
comes slower, i.e, the NEB curve gets flatter, with the 
introduction of boron (Fig. [4|. 

Origin of the reduction in SW activation barrier can be 
explained in terms of changes in bond strength, charge 
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FIG. 5. (color online) Due to the B-doping phono ns shift 
toward lower frequencies, which indicate bond weakening. 
Phonons for the 14-atom pyracyclene around the rotating 
dimer is shown for Ih- symmetry. Arrows are guide to eye. 



density and phonon spectra, which point toward bond 
weakening at and around the rotating dimer. First, we 
calculate C— C, C— B, and B— B bond strengths in pris- 
tine Ceo, C59B, and C58B2, respectively, which can be 
approximately calculated through following expressions, 

E C -c = 

Ec-b = [EcsgB - 87 x E C -c] /3 

£ B -b = [£c 5 sB 2 - 85 x E c -c - 4 x E C -b] , (1) 

where E tot represents the total binding energy of the re- 
spective structures. Calculated bond strength decreases 
with B-doping: 5.87, 4.74 and 3.75 eV, respectively, for 
Ec-c, Ec-b and Eb-b- This hints toward the presence 
of weaker C— B bonds around the rotating dimer, causing 
the SW bond rotation to be easier for the B-doped het- 
erofullerenes. This qualitatively explains the reduction in 
activation barrier. We further note that the amount of 
reduction is related to the reduction in bond strength as 
SE a = n[Ec-c — ^c-b], where n is the number of C— B 
bonds present in C59B (two), and C58B2 (four) hetero- 
fullerenes. 

Next we turn our attention to the phonon frequencies 
calculated at the T point. It is reasonable to assume that 
the atoms around the rotating bond are mostly affected 
due to B-doping and SW rotation. Thus, we considered 
the 14-atom pyracyclene ring around the rotating bond 
to calculate the normal vibrational modes. This signif- 
icantly reduces the size of Hessian matrix (from 180 x 
180 for the full cluster to 42x42) to calculate phonon 
frequencies. Fig. [5] shows that due to B-doping, phonon 
frequencies shift gradually to the lower frequencies, indi- 
cating the bonds get softer as the phonon frequency (v) 
is related to bond stiffness (k) as, v oc \fk. The observed 
softening in phonon frequencies arises due to two factors. 
Firstly, the C— B and B— B bonds are longer and also 
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FIG. 6. (color online) Bonding charge densities Ap for pristine Ceo, and single and double B-doped heterofullerenes. The B 
atoms are shown in green. The Ap for a particular cage is calculated as, Ap = p tot (r) — p° 58 (r) — p x (r) — p Y (r), where X— Y 
represents the rotating C— C, C— B and B— B dimers for Ceo, C59B, and C59B2 structures, respectively. For pristine Ceo, the 
charge accumulation centers (red) appear at the middle of the bond leading to C— C covalent bond, which shifts significantly 
toward the C-atoms (shown with the arrows) due to B-doping leading to ionic character of the C— B bonds. 



less stronger compared to the C— C bonds, as pointed 
out earlier. Additionally, the C— C bonds around the 
rotating unit also get elongated due to the substantial 
strain introduced by B-doping, and thus get weakened. 
The rotation of bond is easier in such softer environment, 
which eventually reduces the SW activation barrier. 

Bonding charge density analysis (Fig. [6| also quali- 
tatively points towards softening of the bonds that are 
involved in SW process. Figure [6] shows the plots of 
bonding charge densities for pristine Ceo, singly and dou- 
bly B-doped Ceo- For pristine Ceo, the C— C bonding is 
completely covalent i.e., the bonding charge accumula- 
tion is at the centers of the bond. This picture deviates 
for B-doped cases, where the charge accumulation centers 
move toward the C-atom for C— B bonds, imparting ionic 
character to the C— B bonds. Thus, the C— C rotation is 
increasingly easier for C59B and C58B2 heterofullerenes. 
The bonding character can be also quantitatively under- 
stood from Bader charge analysis,^ which reveals that 
each B-atom looses ~ l\e\ charge to the neighboring C- 
atoms. This indicate deviation from perfect covalent 
character. 



C. Reaction rate and characteristic time-scale 

At temperatures well below the melting point, the har- 
monic approximation to transition state theory (HTST) 
can be applied to study the reaction rate or character- 
istic time-scale associated with the reaction. The reac- 
tion rate can be expressed in terms of energy and normal 
mode frequencies at the saddle point and initial statef^ 

3iV 

IK 

^HTST = i c -(E a /k B T)^ £ 2 ) 



where E a = (E TS — E l ) is the activation energy with E l 
(E TS ) being the energy corresponding to initial (tran- 
sition) state, and the Vi are the corresponding normal 
mode frequencies. Within the harmonic approximation, 
the the entropic effect to the reaction rate, is included in 
the prefactor involving phonon frequencies calculated at 
zero temperature. The first-order transition state is char- 
acterized by one imaginary phonon, which is excluded for 
the transition state. To calculate the prefactor we have 
only considered the normal modes corresponding to 14- 
atom pyracyclene ring around the rotating dimer. The 
characteristic time-scale associated with a given reaction 
can be calculated from the knowledge of reaction rate 
^htst ag T _ -|y^HTST Figure [7] shows the character- 
istic time required for single ih — » C2V SW transition 
plotted as a function of temperature. The characteris- 
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FIG. 7. (color online) The Ih —> CW transition (reaction) 
rate, and thus the characteristic time scale for the correspond- 
ing transition has Arrhenius dependence on temperature and 
activation energy. Subsequent B-doping in the rotating dimer 
reduces the characteristic time-scale enormously via decrease 
in the activation barrier. 
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tic time scale, measured in second, differs by 10 10 — 10 6 
between Ceo, CeoB and C60B2 in the temperature range 
1000 - 2000 K. For example, at 1700 K, single SW bond 
rotation event takes place in ~ 10 days for pristine Ceo- 
In contrast, at the same temperature, the process is enor- 
mously accelerated and requires only time scale of about 
second and micro-second for C59B and C58B2, respec- 
tively. It should be pointed out that the reverse transi- 
tion (C2V —> ^h, which is not shown in Fig. [7]) is much 
faster than the forward transition (ih — > C*2v, which is 
shown in Fig. [7|. This is caused by the presence of ex- 
ponential factor, exp(— AE/hsT), which is five orders of 
magnitude smaller for C2 V — > ^h, compared to 1^ — » 
at temperature 1700 K. 

IV. CONCLUSIONS 

The Stone- Wales activation barrier in pristine fullerene 
is very large. Earlier it has been found that this high 
barrier can get reduced b y y 3 5% by the presence of an 
extra carbon or hydrogen Our present study shows 
that presence of substitutional boron at the active Stone- 
Wales site can become even more effective in reduction 
of the Stone- Wales activation barrier. It can get reduced 
substantially by ~ 30% and 60%, due to single and dou- 
ble B-doing, respectively. Calculated charge density and 



phonon spectrum indicate that the C— B bonds are softer 
compared to the C— C bonds, and the presence of such 
weaker bonds around the rotating dimer is responsible 
for the observed reduction in activation barrier. From 
thermodynamic point of view, such reduction in the ac- 
tivation barrier AE a enormously reduces the time-scale 
associated with the Stone- Wales process. The reduc- 
tion is of the order of - exp(AE a /k B T) ~ 10 6 (10 12 ) 
fold at temperature 1700 K for single (double) B-doping. 
Stone- Wales defects are known to alter the chemical re- 
activity toward adsorbatespEZI thus are expected to ef- 
fect the proposed reversible H2 storage properties of B- 
het erofullerenes JMlMJ Though, here we have shown the 
reduction of Stone- Wales activation barrier for fullerene, 
similar reduction is expected for carbon nanotube and 
graphene nanostructures due to B-doping, which will cat- 
alyze fusion, and welding process necessary for nanoelec- 
tronic device applications. 
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